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Abstract. In linear elasticity, we decompose the elasticity tensor into two irre- 
ducible pieces with 15 and 6 independent components, respectively. The vanishing 
of the piece with 6 independent components corresponds to the Cauchy relations. 
Thus, for the first time, we recognize the group-theoretical underpinning of the 
Cauchy relations. 

Keywords: linear elasticity, Hooke's law, Cauchy relations 



1. Linear elasticity 

In linear elasticity theory for homogeneous bodies, the stress tensor 
a % i = cr J *, with = 1,2,3, is related to the strain tensor 1 e^i = 

£ik = d( k Uf) by Hooke's law, 

a^ = c^ M e kl , (1) 

see, for instance, [8] or [11]. Here u\ is the displacement field and c l i kl 
the constant 4th rank elasticity tensor. 

Since stress and strain are symmetric tensors, the elasticity tensor 
obeys the symmetries 



c 



ijkl = pjikl _ c ijlk j"2) 



1 We use the notation of Schouten [9] . Symmetrization over two indices is denoted 
by parentheses, A(yj := (Aij + Aji)/2\, antisymmetrization by brackets, := 
(Bij — Bji)/2\. The analogous is valid for more indices, as, e.g., in (13): c 1 ^^ 1 = 
( c »:jfci _ c lk 3 ls j/2\. If one or more indices are exempt from (anti)symmetrization, they 
are enclosed by vertical bars, as, e.g., in (12), c^ k ^ j) = (c lklj + c? kh )/2l, or, more 
complicated, in (6), c (xmm = (c l(kl)j + c j(fci) ')/2 = {c M] + c uk] + c> ku + c? lki )/A. 
Symmetrization over more than 2 indices is, by definition, the normalized sum over 
all possible permutations of the indices involved. Thus, c^"' := (c ljkl + c? %kl + 
c 11 + 21 more terms) /4!, see (5). The totally antisymmetric Levi-Civita symbol is 
denoted by eyjt and e lmn , with £123 = e 123 = +1; for all even (odd) permutations of 
1, 2, 3, we have +1 ( — 1), otherwise 0. The partial derivative d/dx k we denote by dk- 
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A symmetric 2nd rank tensor has 6 independent components. There- 
fore, by collecting the indices i and j into an index pair and k and I 
likewise, the elasticity tensor can be thought of as a 6 x 6 matrix with 
36 independent components. 

Usually one assumes that the stress tensor can be derived from an 
elastic potential W (also called "strain energy function"), i.e., 

a 9W ijkl d 2 W 

OEij OEij OEki 

Then the first and the last pair of indices of c l ^ kl commute, 

c ijkl = c klij ^ 

and the 6x6 matrix is symmetric and carries only 36 — 15 = 21 
independent components. Thus, as is well known, the elasticity tensor 
has 21 independent components. 



2. Generally covariant decomposition of the elasticity tensor 

It is remarkable that the elasticity tensor can be decomposed un- 
der the group of general coordinate transformations (diffeomorphisms) 
into two irreducible pieces. Locally this amounts to the application 
of the 3-dimensional linear real group GL(3,R). Via an analysis of 
Young tableaux, one can determine that a totally symmetric tensor of 
rank r in n dimensions has, expressed as binomial coefficient, ( n+ ^ _1 ) 
independent components, see [9]. Accordingly, 

W c ijkl ._ c (ijkl) /g\ 

has ( 3+ 4 _1 ) = 15 independent components. Because of the symmetries 
(2) and (4), we find 

c (i\(kl)\j) = c i(kl)j = c (i\kl\j) ^ (g) 
and we can go over to the more condensed formula 

{l) c ijkl _ }_(^ c ijkl _|_ c iklj _|_ ^Ijk^ _ ^ 

Alternatively, we can write 

(i) c ijkl = 1 (jjkl + 2c (i\ (kl)\j)^ ^ 

= 1 (jjkl + 2c W^ = 1 (j:> kl + 2c mij) ^j . (9) 
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The surviving irreducible piece of c' J ' fci , with 21 — 15 = 6 independent 
components, encodes the excess of the elasticity tensor over its totally 
symmetric piece, 

(2) c ijkl ._ c ijkl _ (l) c ijkl ^ ^q) 

or, by means of (8), 

(2) c ijkl = 1 ^2c ijM - c iklj - c ilkj ^j = - (c ijkl - c C*IC*=0 U)^ (ii) 

= 1 (cUM _ c i(kl)j^j = 2 ^ijkl _ c (i\kl\j)^j _ ( 12 ) 

We can also express the right hand side in terms of antisymmetric 
terms, 

(2) c ijki = 2 ^ c i[ jk ]i + c i[ji]k^j = 2 (2 c i[jk]i + c i[kl]j^ _ ( 13 ) 

While the tensor c^ 1 ^ has already been introduced in elasticity theory 
[2], this tensor itself is not an irreducible piece of the elasticity tensor. 
In particular, it does not inherit the symmetries of c l i kl . 

The irreducible pieces { a )c^ kl have, as a rule, the same algebraic 
symmetries as c^ kl . We list them here for later use: 

(a) c [ij]kl _ (a) c ij[kl] _ (a) c ijkl _ (a) c klij ^ a = 1,2 . (14) 

Moreover, because of (10) and (5), we have 

(2) c (*jfci) = o. (15) 
Therefore finally we have the decomposition 

c ijkl = (l) c ijkl + (2) c ijkl (16) 

or 

21 = 15 6. (17) 

An analysis by means of Young tableaux guarantees that this decom- 
position is unique and cannot be continued successively. Similar group 
theoretical methods were applied for the decomposition of the electro- 
magnetic constitutive tensor of the vacuum ^ vkX , a 4th rank tensor in 
4 dimensions, see [5]. 
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3. Mapping ( 2 )c'j kl to a symmetric 2nd rank tensor 

Since ( 2 )c y W has 6 independent components, it is to be expected that it 
can be mapped to a symmetric 2nd rank tensor. Clearly, the Levi-Civita 
symbol has to be used in this context. However, the Levi-Civita 
symbol can only be applied to antisymmetric index pairs, otherwise 
information may get lost. Consequently, we should derive, out of ( 2 )c ljfei , 
another tensor which has two pairs of antisymmetric indices. The ob- 
vious method is to antisymmetrize suitably: 

(2) d ijkl ._ (2) c [i\\jk]\l] = (2) c i[jk]l = c i[jk]l _ ( 18 ) 

For the derivation of the last equality, we employed (10). Because of 
(14)i, the antisymmetrization in the definition (18) cannot be applied 
to the first or the last pair of indices. Thus, modulo different conven- 
tions, this is a unique procedure. Note also that the index pairs in (18) 
commute. Thus, we have the symmetries 

(2)-ijkl = _(2)gjki = _(2)gikjl = (2)£lkji = (2) gklij _ 
In contrast to (14)i, 

(2) 5 [<j]fci^ , Wc ij[kl] + 0. (20) 

In order to show that ( 2 )c yfc ' is equivalent to ^c l ^ kl and carries, 
in particular, the same number of independent components, namely 6, 
we have to resolve (18) with respect to ( 2 )c yfc '. We substitute c 1 ^^ 1 = 

(2) d ijkl^ c i[jl]k = (2) d ijlk^ and c i[kl]j = (2) d iklj intQ ^3) and obtain 

(2) c ijkl = 1 (2) d ij(kl) = 2 ^ 2 (2) d ijkl + (2) d iklj^j _ ^ 

This is the inverse of (18) and proves the equivalence of ( 2 )c JJlM and 
(2)£ijkl_ 

We remember that ( 2 )c yM is antisymmetric in il and in jk, see (19). 
Therefore the dual of the "excess" tensor ( 2 )c yfc ' turns out to be 

A-mn := ^rnil^njk ^ ^ . (22) 

The tensor A mn is symmetric and has 6 independent components. 
Indeed, by using (19), we find 

A - L f (2)zijkl _ * (2) ? jilk _ \ (2)zjUk 

^nm ^-nil^rajk <- ^ t mil t njk <- ^-mil^njk c 

= ^mil^njk ^ ^ = A m?1 . (23) 
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With the help of (18), Eq.(22) can be also expressed in terms of the 
original elasticity tensor: 

A mn = ie m iie n jk c ijkl (24) 
= \emiie nj k m c ijkl . (25) 

The last equation is valid, since ^c^ kl is a totally symmetric tensor 
and thus 

~^mil^njk^ = 0. (26) 

In other words, the operator \e m iie n jk extracts the dual of the 2nd 
irreducible piece ( 2 )c^ fc/ from c i: > kl , it is a type of a projection operator. 
Another projection operator is, of course, the symmetrizer (....) in (5). 

Formula (24) has been given earlier in the literature by Haussiihl 
[4]. However, our derivation is new. 

For the inversion of (24), it is very convenient to use the generalized 
Kronecker symbols, see [10]. We multiply the left hand side of (24) by 
two 3-dimensional epsilons: 

pqm rsn a _ }_ pqm rsn ijkl _ }_ <:pq rrs ijkl 

c c i-^mn — ^ c c c mtl c njk — ^ u n u jk 

= * § PQ c i[rs]l = c ]p\[rs]\q] = (2) ff.rsq ^ q e ± ^ 



4. Cauchy relations 

Around 1830, in the early days of modern elasticity theory, Navier, 
Poisson, Cauchy, and others set up molecular models for elastic bodies, 
see Todhunter [12]. If they prescribed, inter alia, certain properties 
to the interaction forces between these molecules that build up the 
macroscopic body, then they found 6 constraining relations between 
the 21 components of the elasticity tensor. Nowadays, they are called 
Cauchy relations. A theory obeying the Cauchy relations and thus 
carrying only 15 independent elastic constants was called rari-constant 
theory, the general case with 21 elastic constants multi- constant theory. 
For decades, there was a scientific battle between both camps. For the 
special case of isotropy, merely 1 independent elastic constant survives 
the Cauchy relations, see Sec. 5. 

More recent discussions of the Cauchy relations can be found, e.g., in 
[1], [3], [4], or [7]. A lattice-theoretical approach to the elastic constants 
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shows, see [6], that the Cauchy relations are valid provided (i) the in- 
teraction forces between the atoms or molecules of a crystal are central 
forces, as, e.g., in rock salt, (ii) each atom or molecule is a center of 
symmetry, and (iii) the interaction forces between the building blocks 
of a crystal can be well approximated by a harmonic potential. In most 
elastic bodies this is not fulfilled at all, see [3]. However, a study of the 
violations of the Cauchy relations yields important information about 
the intermolecular forces of elastic bodies. 

If we take the fundamental meaning of the Cauchy relations for 
granted, then it is near at hand to assume that the vanishing of ^>c % ^ kl 
corresponds to the 6 Cauchy relations: 

(2) c ijki = (Cauchy relations) . (28) 

According to the Young tableau technique, there is no other way to 
split off 6 components from the elasticity tensor in a generally covariant 
and irreducible way. Incidentally, Eq.(28) is, due to (25), equivalent to 
= 0. 

Let us show that (28) really encompasses the Cauchy relations. 
According to (28), (10), and (5), the elasticity tensor is now totally 
symmetric: c l ^ kl = c^ kl \ Consequently, 

c ijki = c ikji ( 29 ) 

In this form, the Cauchy relations can be found in [4]. In "longhand", 
the nontrivial components of Eq.(29) read 

c 1122 =c 1212 j c 1133 =c 3131 > c 2233 = c 2323 j (3Q) 

c 1123 = c 1231 c 2231 = c 2312 c 3312 = c 3123 (31) 

In the 6 x 6 notation, 1 ^ 11 , 2 ^ 22 , 3 ^ 33 , 4 ^ 23 , 5 ^ 31 , 6 ^ 12, 
we find the alternative formulae 

c 12 = c m , c 13 = c 55 , c 23 = c 44 , (32) 

c 14 = c 5 6) c 25 = c 4 6) c 36 = c 45_ (33) 

These coefficients coincide with those of Love in [7] p. 100, Eq.(13). But 
certainly, formula (29) is much more compact and easier to remember. 



5. Isotropy 

Let us test our results for the rather trivial case of isotropy. For an 
isotropic body, the elasticity tensor can be expressed in terms of the 
metric tensor as 

= A gVg" + n (g* k g l 3 + g il g^ , (34) 
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with the Lame moduli A and //, see [8]. Consequently, 

(i) c ijkl = ( A + 2 M ) = + ff <fc^j + gilgj^ . (35) 

The 2nd irreducible piece can be conveniently determined by substi- 
tuting (35) directly into (10): 

(2) c m = ^ gijgk i _ gikg ij _ gUg ^ (36) 

Furthermore, 

(2) d ijki = c i[jk]i = (a - M ) ^b/]' . (37) 
The corresponding 2nd rank tensor (25) turns out to be 

^mn — 1 tmiltnjk 9 9 ~ ^ V • 1<3»J 

Accordingly, 

Ay = — ffij • (39) 

Thus, for an isotropic body, the Cauchy relations translate into A = fi. 

Then the elasticity tensor depends only on one modulus. For most 
isotropic elastic bodies this is certainly not fulfilled. 



Appendix 

An analogous decomposition can be made for a 4th rank tensor d tjkl 
which lacks the "Onsager symmetry" (4) but which still carries the 
symmetries (2) 

d ijki = d jiki = d ijik _ ( 4Q ) 

The piezooptical tensor is of such a type, see [4] p. 269, for instance. 
The first irreducible piece is defined as in (5): 

Wjijkl ._ d (ijkl) _ 1 ( jijkl _j_ jikjl _|_ jiljk 

6 V 

+d kHj + d 3Uk + jjkilj _ ( 41) 

Then we quantify the violation of the pair symmetry by 

(3) d ijki ._ 1 ^pjkl _ d kiij^ ( 42 ) 

and the rest as 

(2) d ijkl ._ d ijM _ (l) d ijkl _ (3) d ijkl (43) 
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Finally we have the decomposition 

^ijkl _ {l)^ijkl _|_ (2)^ijfcZ _|_ (3)^ijM 

or 

36 = 15 ©6© 15, (45) 

which is irreducible. 

In analogy to (22), one can map ( 2 )ff-? fc ' to a 2nd rank tensor: 

A mri := -^mil^njk ^ ^ d^ . (46) 

It is again symmetric, that is, A[ mn j = 0, as can be shown in an 
analogous way as in (23). Moreover, since = ^d^ kl \ we have 

\tmiitnjk ^d ljkl = 0. As for the third part, 

~^mil^njk^ ^ d ^ = —£mli^nkj ^ ^d ^ = —^mil^njk^ ^d "' 

= -\e mil e njk ^d^ kl = 0. (47) 

Consequently, 

A mri = -^rnil^njk d ^ ■ (48) 

Imposing the additional symmetry d l ^ kl = d khj , that is, ^>d^ kl = 0, 
leads back to the results found for the elasticity tensor c^ kl . 
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